Exclusive production of quarkonia as a probe of the GPD E for gluons 
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Exclusive quarkonium photo- and electro-production off the nucieon is studied in the framework 
of generalized parton distributions (GPDs). The short distance part of the process is treated at 
leading order in perturbative Quantum Chromodynamics. The main focus is on the GPD E for 
gluons. On the basis of different models for we estimate the transverse target spin asymmetry 
for typical kinematics of a future Electron Ion Collider. We also explore the potential of measuring 
the polarization of the recoil nucieon. 
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I. Introduction. — For about 15 years, GPDs [H-Q have 
been playing a key role in hadronic physics for a number 
of reasons. First, GPDs serve as unifying objects, con- 
taining the information encoded both in ordinary parton 
distributions and in form factors. Second, GPDs allow 
us to explore t he p artonic structure of hadrons in three 
dimensions 13" 3|- Third, GPDs enter Ji's spin sum rule 
of the nucieon 3| ■ 

GPDs can be measured in hard exclusive processes like 
deep- virtual Compton scattering off the nucieon or hard 



They depend 



exclusive meson production [2t5, 
on three kinematical variables: the average momentum 
fraction x of the partons, the longitudinal momentum 
transfer ^ to the nucieon (skewness), and the invariant 
momentum transfer t of the process, i.e., F = F(x,S.,t) 
for a generic GPD F. According to [2|, GPDs give access 
to the angular momenta of quarks and gluons inside the 
nucieon, where the total spin of the nucieon is given by 
^ = J'' + . Specifically, the gluon angular momen- 
tum can be determined through ^ 



1 

- dx(^Haix,C,0) + E\x,^,0)), (1 



with and i?^ denoting the dominant (leading twist) 
GPDs of unpolarized gluons inside the nucieon. Consid- 
erable information on is already available, for it is 
connected to the ordinary unpolarized gluon distribution 
via H^{x, 0, 0) — xg{x) — see for instance Refs. Til- 22|. 
In comparison, our knowledge about is still marginal. 
Therefore, in particular, the value for in Eq. ([T]) is still 
very uncertain. 

It has been known for quite some time that exclusive 
quarkonium ( J/V" or T) production off the nucieon is very 
suitable for probin g th e gluonic structure of the nucieon 
in a clean way (23l. 124^ . since quark exchange plays only 
a minor role. Moreover, due to the large scale provided 
by the heavy quark/meson mass, perturbative Quantum 
Chromodynamics (QCD) can be applied even for photo- 
production. In the present work, we consider both photo- 
and electro-production of J/i) and T off a proton target. 



Using leading order (LO) results for the hard scattering 
coefficients we study the prospects for measuring E^ by 
means of quarkonium production. To this end we con- 
sider several models for , and compute the transverse 
target spin asymmetry as well as a double spin observ- 
able which requires polarimetry of the recoil nucieon. We 
provide numerical results for typical kinematics of a po- 
tential future Electron Ion Colfider (EIC) [25l-[27|. 
II. Theoretical framework. — For definiteness, we con- 
sider the process 



7^*' (9, li) + p(p, v) ^ V{q', //) + p{p', v'), 



(2) 



where the 4-momenta and the helicities of the particles 
are specified. We further use = —q^, vn? = = p'"^ , 
my — q''^ , t — [p ~ p'Y , and the squared photon- nucieon 
cm energy = (p + qY ■ The skewness variable can be 
expressed as 



Xb 



Xb 



with Xb 



(3) 



which holds for arbitrary values of Q^. The minimal 
value of t is given by |io| — 4to^^^/(1 — ^^). 

For much larger than all other scales of the pro- 
cess, an all order proof of QCD factorization has been 
formulated for the process in ([2]) [l3|- In the case of 
quarkonium production, one may expect factorization to 
hold for arbitrary Q^. In fact, a next-to- leading order 
(NLO) calculation of the unpolarized photo-production 
cross section is compatible with factorization (isj . 

To LO in the strong coupling there are six Feynman 
graphs — see Fig.[T]for a sample diagram. They factorize 
into the hard subprocess 7*3 — >■ Vg and gluon GPDs. We 
computed the helicity amplitudes of the subprocess in the 
non-rclativistic approximation for which the heavy quark 
and antiquark carry the same momentum. Our results 
agree with previous calculations 28-33- The structure 
of the LO amplitudes im plie s that one is sensitive only 
to the GPDs H3 and 0,113 . 

For transversely polarized photons and vector mesons. 
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FIG. 1: Sample LO diagram for the process in ([2]). The lower 
part of the diagram is parameterized in terms of gluon GPDs. 



the nonzero helicity amplitudes A^^v'./ji/ of the full pro- 
cess read 




(4) 
(5) 



{x + ^){x - ^ + ie) 



for a generic GPD F. In Eqs. (U), ([S]) we use the def- 
initions H^g = H9 - ^2/(1 - (^)E3, t' = t - to, and 
C ~ Wneqeasfvmv / {Nc{Q^ + my)), where fv denotes 
the quarkonium decay constant. Moreover, one has 



Q 

Mov'fiu = A^±i/',±y 

mv 



(6) 



for longitudinal transitions. A corresponding relation 
between the longitudinal and the transverse amplitudes 
was previously obtained in the pioneering work on exclu- 
sive J/^ production in the leading double- log approxi- 
mation [23|. 

III. Generalized Parton Distributions. — For the GPD H 
we take t he p arameterizations obtained in previous anal- 
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that quark GPDs are also needed 
since the GPDs are evolved to different scales. In the 
case of E, we use the valence quark distributions from 
Ref . [Sll [S^l , while we explore different scenarios for glu- 
ons and seaquarks. They are modelled through double 
distributions [l|, [s^ according to 



E\x,(„t) 



d(3 



1-1/31 



f\P,a,t) = E\P,Q,t) 



1 + 1/31 

15 [(1 



da5{P + ^a-x)P{l3,a,t), 



a 



212 



16 (i-mf ' 

with £:'(^,0,t) =e^=*|;3|-<*£;*(|;3|,0,0). 



(7) 



For gluons we define xe^{x) = E^{x, 0, 0) and investigate 
the two cases 



e^{x) = N^x-'^-^' 
e^ix) = N^x-^-^' 



(1 — tanh(l — x/xo), 



(8) 
(9) 



where the ansatz in Q has a node at a; = xq. Such 
a possibility is currently not ruled out We further 



Var. 




ATS 


Xo 




J" 


N" 


r 


1 


0.15 










0.214 


-0.009 


0.014 


2 


0.15 


-0.878 




-0.164 


0.132 


0.156 


0.041 


3 


0.10 


-1.017 




-0.190 


0.119 


0.182 


0.045 


4 


0.10 


3.015 


0.05 


-0.190 


0.119 


0.182 


0.045 


5 


0.10 


-1.974 


0.3 


-0.190 


0.119 


0.182 


0.045 



TABLE L Parameters of and e« at the scale ^ = 2 GeV. 
For gluons, the Variants 1,2,3 refer to the ansatz in ((8|, while 
Variants 4,5 refer to ®. Also shown is the second moment 
620, and values for the angular momenta as defined in U]). 



Ei{x,0,0) 





, and use a flavor-symmetric 
= e^. For e^{x) we make 



define e'^{x) 
sea, i.e., ; 

an ansatz analogous to ([8|) and do not consider a node. 
For the parameters be, a'^, and 5e we do not distinguish 
between gluons and seaquarks. 

Two constraints have to be satisfied when fixing the 
parameters for and . First, the momentum sum rule 
for unpolarized parton distributions in combination with 
Ji's spin sum rule leads to 
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^20"' 2 ^ e| 



20' 



with 



dxx'^-^e^x). 



(10) 
(11) 



Second, the density interpretation of GPDs in the im- 
pact parameter space [isl leads to a positivity bound for 
and e^ — see Refs. [19, 3Q-32] for more details. In 
particular, one finds b^ < bh, a'^ ^ oi'h^ where bh and a'^ 
appear in the double distribution ansatz of . We take 

-J- (with /X repre- 



0.25 GeV~^ In ■ 



2.58 GeV" 

senting the scale of the GPD) and a'y^ — 0.15 from pre- 



vious work l3,l20j. We choose bf. — 0.9 bh and explore 
two different values for a'^ (see Tab. [J). Moreover, we use 
6" = 0.1, as well as ;3f = 6 and /3| = 7 H. After these 
choices have been made, only the normalization constants 
and remain to be determined. 
We parameterize and at the scale /i = 2 GeV. For 
we consider five different variants, where the respective 
parameters are listed in Tab. [J Variant 1 means = 0, 
and the normalization N'' is fixed by means of the rela- 
tion in (jlOp . (There is actually some support for a rather 
small E^: this GPD has a model-dependent relation to 
the transverse momentum dependent gluon Sivers func- 
tion [3^, which may be small (36l - l38| .) In the remaining 
four cases we first determine by saturating the pos- 
itivity bound, then compute from (fTO|) . and finally 
check whether satisfies the positivity bound. Variants 
4 and 5 for contain a node. The positivity bound does 
not allow one to fix the sign of N^. However, none of our 
general conclusions depends on this sign [30|. We also 
checked that all variants for and are compatible 
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;x = 3.1GeV, ^ = 5 x 10-^ t = ^ = 9.46GeV, 5 = 5 X 10-^ t = 
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FIG. 2: Variants 1-5 for , together with H^, at the scales 
Pl = mjj^ = 3.1 GeV (left) and jU — mr — 9.46 GeV (right). 



FIG. 4: SSA An in ^ for photo-pro(iuction of J/t/) (left) 
and T (right) as function of W for different variants of . 
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FIG. 3: Density in gU for Variant 2 (left) and Variant 4 
(right) at a; = 0.005 and jj, — 2 GeV. The outer ring indicates 
half of the maximum density. 



with a preliminary data point for the transverse target 
spin asymmetry for exclusive (f) production from HER- 
MES [Sflll^. When calculating observables we evolve 
the GPDs to the scale fi — {my + Q^)^/^ by using the 
code of Vinnikov ji^l- In Fig. [21 the GPDs are displayed 
at the scale of the quarkonium masses. 

Following 'i'l, the gluon and strange quark contribu- 
tions to the nucleon spin are = (/ijo + ^20) 1"^ '^^^^ 
as J* = (/120 + 620), where the densities /i^/^ are related 
to iJf /« in the same way as e^/* are related to Tak- 

leads to 



ing /ifo = 0.4276 and /i|o = 0.0153 from [32,1411 
the values for and J" in Tab.[ll Because of (jlOp . a 
change of J* implies a change of . For our paramc- 
terizations, the contribution from to the nucleon spin 
can be significant (up to about 20 %). 

One can also calculate the density of unpolarized glu- 
ons in transverse position (impact parameter space. 
If the nucleon is transversely polarized (along the X- 
direction), this density is given by (l3| 



bl_d_ 

m dbl 



£%x,bl), (12) 



with W and £^ denoting the &_L-space representation 
of and , respectively. The sample plots in Fig. [31 



show, in particular, how much the maximum of the den- 
sity in is shifted away from the origin due to the 
presence of the i?^-term. 

IV. Polarization observables. — As discussed in Sect. II, 
one has two independent amplitudes: the non-fiip am- 
plitude A^++,++, dominated by H^, and the spin-flip 

amplitude Mj, which is determined by E^ . The 

following observables allow one to measure those ampli- 
tudes (up to an overall phase): the unpolarized cross sec- 
tion, the transverse target single spin asymmetry (SSA) 
An (with the polarization being normal to the reaction 
plane — often is also denoted as Aut), and two dou- 
ble spin asymmetries (DSAs) requiring a polarized target 
and polarimetry of the recoil nucleon ^3ft] . Here we focus 
on 



An 
Als 



-Ml 



-) 



\M++,++\^ + 
2Re(7W++,+^ 



\M+- 
M 



+-.++' 



\M++,++\^ + \M+-^++\ 



(13) 
(14) 



where for the DSA Als the target is longitudinally polar- 
ized, and the outgoing nucleon is transversely polarized 
(in the reaction plane, "sideways") [sO]. We consider pro- 
duction of both J/ip and T for typical EIC kinematics. 
While the J/ip final state has a larger cross section, one 
can expect a better convergence of the CKs-expansion in 
the case of the T [i^ . A detailed comparison with exist- 
ing data for the unpolarized cross section will be given 
elsewhere fsO*] — see also Ref. [28j . 

In Fig.m An is shown for photo- production of J/ip and 
T as a function of W. This asymmetry is rather small for 
most variants of E^ , mainly because the respective non- 
flip amplitude and the spin-flip amplitude have a similar 
phase. It can also be seen, however, that larger values 
of An are currently not ruled out. On the basis of a LO 
calculation one can not draw a definite conclusion about 
the optimal W for a measurement of An- But higher 
order terms to the unpolarized cross section are better 
under control for lower values of W [2^. In general, the 
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FIG. 5: DSA Alt in (HU for photo-production of J/i/- (left) 
and T (right) as function of W for different variants of . 
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FIG. 6: An (left) and ^ls (right) for production of J/il) as 
function of for different variants of . 



spin asymmetries are less influenced by NLO corrections 
than the cross section The DSA Als is displayed 

in Fig. [5] This observable is small only if is small. It 
is worthwhile to explore the feasibility of a correspond- 
ing measurement, since A^s may give the most direct 
access to E^. (We note that, from a theoretical point of 
view, the DSA Asl is equally well suited js^].) Finally, 
the Q^-dependence of A^ and Als for J/ip production 
is shown in Fig.lHl Electro-production at low is at- 
tractive because of the large count rates. For > my, 
higher order corrections may be less important [i^j, but 
an explicit calculation does not yet exist. 
V. Summary — We have explored the potential of mea- 
suring the GPD E^ through exclusive production of 
quarkonium at a future Electron Ion Collider. The study 
is based on a LO calculation of the short distance part 
of the process, and several models for E^ which respect 
the currently known constraints. Most variants of E^ 
lead to a rather small transverse target SSA Ajy, but a 
healthy An is presently not ruled out either. We have 
also found promising results for a double polarization ob- 
servable (polarized target and polarized recoil nucleon), 
which provides a quite direct access to E^ . 
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